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The exclusive charge exchange reaction pD —> n(pp) at intermediate and high ener¬ 
gies is studied within the Bethe-Salpeter formalism. The final state interaction in the 
detected pp pair at nearly zero excitation energy is described by the 1 S'o component 
of the Bethe-Salpeter amplitude. Results of numerical calculations of polarization 
observables and differential cross-section persuade that, as in the non-relativistic 
case, this reaction (i) can be utilized as a “relativistic deuteron polarimeter” and (ii) 
delivers further information about the elementary nucleon-nucleon charge-exchange 
amplitude. 


I. INTRODUCTION 

The investigation of polarization observables in electromagnetic and hadronic processes 
at high energies provides refinement of the information about strong interaction at short 
distances and the relevant reaction mechanisms. Accordingly, the experimental study 
of processes with polarized particles becomes more and more important. Experiments 
with deuteron targets or beams B0BB are particularly interesting, since the deuteron 
serves as a unique source of information on neutron properties at high transferred mo¬ 
menta; the knowledge of which allows, e.g. to check a number of QCD predictions and 
sum rules. For example, for an investigation of the NN interaction in the deuteron at 
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short distances, the three deuteron form factors (magnetic, electric and quadrupole) have 
to be determined. In the elastic eD scattering with unpolarized particles one can measure 
only two independent quantities, e.g. the magnetic form factor and the deuteron func¬ 
tion A(Q 2 ), the latter being a kinematical combination of all three form factors. Even 
these two quantities reveal an important information about the quark physics and dy¬ 
namics at short distances as demonstrated, for instance, in recent measurements [1] at 
TJNAF. However, for a full determination of the deuteron form factors separately, one 
needs measurements with polarized particles. For instance, measurements of the tensor 
analyzing power T 2 o of recoil deuterons in elastic eD scattering allow for a determination 
of the charge form factor G c at high transferred momenta. Namely the charge form fac¬ 
tor is very sensitive to details of the NN interaction 0. [3], |5|] and, besides information 
about short range correlations in the deuteron, the investigation of G c may essentially 
constrain the theoretical models applied in this area. However, in spite of the fact that 
the electromagnetic processes are considered as the cleanest ones, in such reactions one 
probes mainly the quark structure of the target, leaving almost untouched the physics 
connected with gluon degrees of freedom. In this context, hadron deuteron processes can 
be considered as complementary tool in investigating phenomena at short distances and 
also as a source of unique information unavailable in electromagnetic reactions, such as 
nucleon resonances, checking non-relativistic effective models, NN potentials etc. 

Experimental and theoretical investigations of the proton-deuteron processes at inter¬ 
mediate and high energies have started some decades ago by studying elastic pD scat¬ 
tering ||, exclusive and inclusive break-up reactions [TL |3|] with the goal of determining 
the the details of the deuteron wave function at short distances and the relevant reac¬ 
tion mechanisms. Note, that in elastic backward pD processes it is possible to determine 
completely the reaction amplitude by measuring a full set of polarization observables (see, 
e.g„ gEJEl). Hence, as in the electromagnetic processes one needs to measure different 
polarizations of the recoil deuteron. 

Since such polarization observables can be studied only by an additional secondary 
scattering of the reaction products, e.g., inside a polarimeter, it is obvious that second 
process must possess a high enough cross section to assure a good efficiency of the po¬ 
larimeter. Traditionally, at low energies one uses the process 3 He (D,p) 4 He Jl2|, while 


for relativistic energies the elastic Dp scattering serves as polarimeter [13]. Bugg and 
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Wilkin [[Uj proposed to use as an efficient deuteron polarimeter the process pD —> n(pp), 
where the hnal pp pair is detected with very low excitation energy (see also ref. !§)• it 
has been argued that at low excitation energies and low transferred momenta the process 
pD —> n(pp) is determined by the elementary pri charge-exchange of the incoming proton 
with the neutron in the deuteron, whereas the second proton acts merely as a spectator. 
In this case, the detected pp pair can be considered to be in the 1 S'o final state. Within 
the non-relativistic spectator mechanism, a significant value of the tensor analyzing power 
T 2 o and a vanishing vector analyzing power were predicted [|F|] . The corresponding cross 
section is rather large, so that the process pD —> n(pp) can be considered as a good 
tool for determining the deuteron tensor characteristics. Later detailed investigations 
[ |l6l , [I?]- 0, [HI of this process confirmed the previous theoretical predictions |L4 


The 


charge-exchange processes of this type have also been proposed for investigations of other 
processes with deuterons, e.g. pp —> Dn + reactions [0, AN systems [0, NNtt systems, 
inelastic (D, D') reactions off heavy nuclei to study isoscalar transitions AT = 0, AS = 1 
etc. In addition, the reaction pD —> n(pp) may be interesting in investigations of the 


elementary NN amplitude. As shown in ref. 0 , this reaction can be used as a part of 
a complete set of experiments to determine completely the amplitude of charge-exchange 
reaction pn —» np fl23|. One may also investigate the influence of the nuclear medium 


on the NN amplitude P3L E3, or to study the double spin flip processes in quasi elastic 


scattering of deuterons from heavy nuclei |I6| 


The direct consequence of these facts is that nowadays the interest in investigations 
of charge-exchange processes does not abate. For example, at the cooler synchrotron 
COSY in FZ Jiilich a program to study of similar processes at relativistic energies has 
already started [0, ||] and a detailed investigation of polarization observables is envisaged. 
Inspired by this, in previous work |26|] we investigated the process pD —> n(pp) within 
the impulse approximation. The goal of the present paper is to consider the charge- 
exchange reaction at relativistic energies, as accessible at COSY and upgraded Dubna 
accelerator by (i) taking into account the effects of hnal state interaction, (ii) to check 
whether in this case the non-relativistic predictions 13 hold and (iii) whether the reaction 
still can be regarded as a deuteron polarimeter tool and/or as complementary source of 
an experimental determination of the elementary charge-exchange partial amplitudes. 
We propose a covariant generalization of the spectator mechanism 0 based on the 
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Bethe-Salpeter (BS) formalism and on a numerical solution of the BS equation with a 
realistic one-boson exchange kernel |27], |2^| . Our amplitude of the process in an explicitly 
covariant form allows for a determination of any polarization observables. Nevertheless, 
here we focus on calculations of the cross section and the tensor analyzing power T 20 for 
kinematical conditions as relevant for experiments at COSY. Vector analyzing powers of 
the deuteron are strictly equal to zero in our case, since we consider the 1 S'o NN final 
state. The corresponding expressions are quite lengthy, so we do not present an explicit 
comparison with the non-relativistic formulae, nevertheless, the final results are written 
in a form as close as possible to the non-relativistic case. We are going to compare our 
results, computed at non-relativistic energies, with corresponding data and corresponding 


non-relativistic calculations ||T8| thus demonstrating that our formulae hold in the non- 
relativistic limit as well, as it should be. We adopt for the process pD —> n(pp) with 
a slowly moving, small-excitation energy pp pair the same mechanism as in [|14|], i.e., 
the process is treated, in the deuteron center of mass, as a charge-exchange between 
the incoming proton and internal neutron with the second proton as a spectator. The 
resulting pp pair is supposed to be detected solely in the 1 S , 0 state. Particular attention 
is paid to the pure impulse approximation, where the final state interaction in the pp 
pair is, for the time being, disregarded. Within the impulse approximation we study 
systematically peculiarities of the reaction to find proper kinematical conditions where 
the supposed mechanism is adequate and to fix the choice of input parameters. Then the 
final state interaction in the 1 Sq state is taken into account by solving the inhomogeneous 
BS equation in the one-iteration approximation, which allows one to numerically compute 
the corresponding partial amplitudes. The effects of final state interaction are found to be 
substantial and essentially improve the agreement with data. Methodological prediction 
for the COSY kinematics will be presented as well. 

This paper is organized as follow. In Sect. 11 kinematics and notation are introduced. 
A short review of the spin structure of the amplitude — + 1 —> — + 0 and definitions 
of polarization observables are presented. Sect. Ill deals with the invariant amplitude 
within the BS formalism and defines the corresponding BS wave functions. A reduction 
of the covariant form of the amplitude to the traditional form in the two dimensional 
spinor space is also performed in this section. In the next section IV a detailed study 
of relativistic impulse approximation and a comparison with data is given. In Sect. V 
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the procedure of accounting for final state interaction in the continuum within the BS 
formalism is discussed. The total BS wave function for the 1 So configuration is presented, 
the corresponding numerical calculations of the cross section and tensor analyzing power 
and a comparison with the available experimental data is performed. Conclusions and 
summary may be found in Sect. VI. Some cumbersome expressions are relegated to the 
Appendices. 


II. KINEMATICS AND NOTATION 

Since in the exclusive processes with three nucleons in the final state we are interested 
in studying correlations in the pp pair we select those of them which, in the deuteron 
center of mass system, correspond to final states with one fast neutron and a slowly 
moving proton-proton pair, i.e. reactions of the type 

p + D = n + (pi + p 2 ). (1) 

A peculiarity of the processes (HD is that the transferred momentum from the proton to the 
neutron is low, hence the main mechanism of the reaction can be described as a charge- 
exchange process of the incoming proton off the internal neutron, whereas the second 
proton in the deuteron remains merely as a spectator. As well known (see e.g. ref. f29|| ) 
the differential cross section of elementary charge-exchange process pn —> np exhibits a 
sharp maximum at vanishing transferred momenta. Therefore, if the reaction (jl|) is indeed 
governed by a charge-exchange subprocess, then the resulting pp pair will be detected with 
low total and relative momenta. Reactions of this kind can fairly well be distinguished 
from other processes. For relatively low initial energies, such reactions are quite well 
experimentally investigated. In Fig. [lj. the diagram of such processes is schematically 
depicted. The following notations are adopted: p = (E p , p) and n = (E n , n) are the 
4-momenta of the incoming proton and outgoing neutron, P' is the total 4-momentum 
of the pp pair, which is a sum of the corresponding 4-momenta of detected protons, 
Pi = (Ei, p^, P2 = (E 2 ,p 2 ), P' = Pi + p 2 . The invariant mass squared of the pair 
is Sf, Sf = P ' 2 = (2m + E x ) 2 , where m stands for the nucleon mass and E x for the 
excitation energy of the pair. Conform the supposed reaction mechanism, the excitation 
energy E x ranges to a few MeV, say E x = 0 — 8 MeV. At such low values of E x , the 
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main contribution in the final state of the pp pair in the continuum comes from the 1 S'o 


configuration |L8| . In what follows all corrections from higher partial waves are neglected, 
however, we realize that for higher values of E x an increasing role of these corrections is 
expected. 

Further, the Dirac spinors 


normalized as u(p)u(p) 
reaction m reads 


/ 


u 


(p,r) = \Jm + 


Xr 

> ' P) 


( 2 ) 
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V m + e 

2 77i are introduced. Then the differential cross section for the 
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where A (p, D) is the flux factor, Mjj is the invariant amplitude and the statistical factor 
1/2 is due to two identical particles (protons) in the final state. By changing in eq. (|3|) 
the kinematical variables from the momenta 2 to the relative and total momenta of the 
pair and taking into account that there is no angular dependence in the 1 Sq configuration 
the cross section can be written as 


d 3 cr 


Am 2 


\M„\ 


d 3 n 


16tt^/A0 o,D) \ s f 1 2P n (2vr) 3 2 


( 4 ) 


Since our numerical solution of the Bethe-Salpeter equation has been obtained in the 
deuteron rest system all further calculations are performed in that. The quantization z 
axis is chosen along the momentum p of the incoming proton; the x and y axes will be 
specified below. Changing the variables in (j3|) we arrive at 


d 2 a 
dt dsj 


1 1 
2 647rA(p, D) \ 


1 - 


4 m 2 
s f 


d(j> 

(2tt) 3 


\M,P 


( 5 ) 


where q = n — p, t = q 2 , Sf = (D — q) 2 , and 0 denotes the azimuthal angle of the final 
neutron. Further we consider only the case where the initial proton and the final neutron 
are unpolarized and the polarization density matrix of the initial deuteron, p£>, possesses 
an axial symmetry relative to the ^-direction, i.e. 

Pd = ^ 1 + Pv Tio + Pt T 2 o, 
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where p v and p t are the vector and tensor polarization parameters, respectively. In this 
case the angular dependence upon 0 in eq. (Q) is trivial. Finally one has 


d 2 a 
dt dsj 


1 1 
2 647rA(p, D) \ 


1 - 


4m 2 1 

s f (2vr) 2 


\M fi \ 2 


( 6 ) 


where the amplitude \Mfi\ 2 can be computed at arbitrarily fixed value of 0, e.g., 0 = 0. 
The procedure of computing the amplitude Mfi consists on several stages: (i) we analyze 
general spin structure in terms of a decomposition of Mfi over a relevant independent set 
of spin variables with coefficients being invariant partial spin amplitudes of the process 
since all observables can be expressed via these partial amplitudes, (ii) the diagram in 
Fig. I is computed explicitly and the results for Mfi are regrouped to obtain an expression 
in the same form as for the general decomposition of Mfi , (iii) from the direct comparison 
of the expression with the phenomenological form the partial spin amplitudes are found 
and the polarization observables computed. 

By virtue of zero angular momentum of the final pair, the process ([!]) is of the type 
1/2 + 1 = 1/2 + 0, for which the symmetry restrictions leave only six independent 
(complex) partial amplitudes. The choice of their explicit representation depends upon 
the kinematical conditions of the attacked problem. One may choose the hclicity repre¬ 
sentation, or the representation with a given spin projection for specific choices of the 
quantization axis etc. In this paper, we choose the following way to determine the par¬ 
tial amplitudes (see also [0): initial f) and final |/) states of the system, besides other 
quantum numbers, are characterized by the spin projections on the z axis; in the matrix 
element Mfi this spin dependence is written explicitly by emphasizing in |i) and |/) the 
3-polarization vector of the deuteron and the two-component Pauli spinors for nucleons. 
We introduce, in the deuteron center of mass, where D = 0, the three basis vectors as 
follows 


p_ Ipxnj 

|p|’ I [p x n] 


a = [b x c]. 


Then the amplitude Mfi can be represented in the form 


(7) 


Mfi = TX = [xt']a (Map £m) [Xr]p, «,/?=!, 2, 


( 8 ) 


where r',r and M are the spin projections for the neutron, proton and deuteron, respec¬ 
tively. The amplitude M a p is a vector in the coordinate space and a matrix in the spinor 
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basis and consequently can be decomposed over the introduced basis vectors (0) and Pauli 
matrices cy (i = x, y, z) as 


M a p — i A b 5 a p + B b (<r • b ) a p + Ca(cr • a) Q/ g + 

+ V a (<r • c) a/3 + 5 c (<r • a) Q/3 + T c (cr • c) a/3 . (9) 


In eq. ([ 

system: 


$, M stands for the polarization vector of the deuteron in its center of mass 


1 

f 1 ) 

1 

( i ^ 
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( 10 ) 


Due to the use of the two-dimensional spin and 3-dimensional vector representation, eqs. 
(||) and (D are not manifestly covariant. Nevertheless, such a representation of par¬ 
tial amplitudes is of most general form and valid in both relativistic and non-relativistic 
considerations. This may immediately be seen if one expresses in covariant matrix ele¬ 
ments the polarization 4-vector of the deuteron (in any system of reference) via the 
3-dimensional as 




(D ' £m) + , (D ' £m) 

Mr, ’ M Md(Ed + Mb) 


and passes from the 4-spinors defined in eq. (Q) to Pauli spinors Xr- 

For convenience, the y and x axes are oriented along b and a, respectively. In this 
case the neutron azimuthal angle (j) can be put equal to zero. The z axis, as mentioned 
above, is parallel to c. The three variables, upon which the partial amplitudes A, B, ...T 
depend, are chosen to be the total initial energy s, the transferred momentum t, and the 
invariant mass of the pp pair Sf. These amplitudes are related with the spin amplitudes 
TX via the following expressions 


A = ( T \_! + T 1 \)/ V 2 , B = - T) 1 ! i)/a/2, 

2222 22 22 

C = -XLi + T. 1 ! i)/V2, v = (T_ 1 i _i - Th )/ V2 , 


£ = T,_ i 

2 2 


T = T[\ 
2 2 


( 11 ) 


Note that having computed the amplitudes A, B, ...JF, all polarization observables for the 
process ([[]) can be found as proper combinations of these partial amplitude. So, if an 
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operator O corresponds to a measurable physical quantity then its mean value is 

^ ^ 6 Ty(MM+) ’ ^ ^ 

where the denominator corresponds to the cross section of the reaction ([!]) with unpolar¬ 
ized particles 


: Tr (MM~ 


= -(AA* + BB* + CC* + VV* + ££* + TV* 

O 


For instance, the tensor analyzing power (T 20 ) is 

Tr (MT 2Q M + ) 


{T20} — 6 - 


Tr (M M+ 

V2 


(AA* + BB* + CC* + VV* - 2 [££* + TT*\). 


(13) 


Tr (M M- 

Note that the representation of the amplitude Mfi by eqs. (§) and @ holds if initial and 
final states can be described by wave functions (pure spin states), otherwise for mixed 
states the square of Mfi must be averaged with the spin density matrices 


= 5X»H P D = l -Y 1 \M){M\. 


(14) 


M 


An explicit covariant expression for the deuteron density matrix can be found in ref. 
[[01] . In the present paper the amplitudes Tjf and observables for the reaction ([I]) are 
evaluated within the BS formalism. It is worth noting that in theoretical considerations 
of high energy reactions with deuterons as target and two interacting nucleons in the final 
state (scattering or bound state) always, at least one of two-nucleon systems is moving, 
consequently Lorenz boost effects must be treated in a consistent way. In our opinion, 
the most appropriate approach for these purpose is the BS formalism where a consistent 
description of the deuteron bound state and scattering states of the NN pair as well as the 
off-mass shcllness of nucleons and Lorenz boost effects may be achieved |[TT] . There are 


other approaches to the relativistic description of reactions with deuterons. For instance, 


the Gross equation [|31|, which is a variant of the BS approach with one nucleon on-mass 
shell; it provides a covariant description of processes like ([[]). An analysis of results within 
the Gross and BS approaches shows that for internal relative momenta up to |k| ~ 1.5 
GeV/c the deuteron amplitudes and wave functions are almost identical. 
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In what follows we compute the amplitude by evaluating the diagram in Fig. |l| 
within the BS formalism. Neglecting the initial state interaction between the incoming 
proton and the deuteron and the final state interaction of the outgoing neutron with the 
pp pair simplifies the calculation. The initial and final states can then be written as 
direct products of spinors of the fast particle and the BS amplitudes of the NN system, 
which correspond to solutions of the BS equation for bound (i.e. the initial deuteron) or 
scattering states (i.e. the pair in the continuum). 


III. INVARIANT AMPLITUDE 


By using the Mandelstam technique |33[ the covariant matrix element corresponding 
to the diagram in Fig. [L] can be written in the form 


D 


TX = uZ ( n ) u r 5 (p) / d 4 k $ P ,(1 + k) af} (— + k — m) a ^ M (k), u AZ Sl/ . 


(15) 


In eq. (|I5D the deuteron BS amplitude $m and the conjugate amplitude <hp/ of the pp 
pair are solutions of the corresponding BS equation, and the charge-exchange vertex A ce 
corresponds to a 4-point Green function of the subprocess pn —> np with, in the most 
general case, off-mass shell nucleons. 

It is convenient to change from outer products of spinors and amplitudes to the usual 
matrix structures. For this sake we redefine the BS amplitude f34| as 


${k) = *{k)U c , #(fc) = 7o^)7o, 

where Uc = ^7270 is the charge conjugation matrix. Then the new amplitudes T (k) may 
be considered as usual 4x4 matrices in the spinor space, and the BS equation becomes 
an integral matrix equation. To find a numerical solution of the BS equation usually 
the amplitude is decomposed over a complete set of matrices, and one solves the resulting 
integral equation for the coefficients of such a decomposition. These coefficients are known 
as the partial BS amplitudes. There are eight independent partial amplitudes for the 
deuteron, and the specific form of them depends on the chosen matrix representation. 
In the present paper we choose the p-spin representation [[k| for the partial amplitudes 
and, since in the considered reaction the transferred momenta is rather small, all the 
amplitudes with at least one negative p-spin are disregarded as they may play a role 
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only at high momenta (see, e.g. [|36[ for the justification). Then we are left with two 
“++” partial amplitudes known as S and D components within the p-spin classification. 
The numerical solution for the BS amplitudes has been found ^7\, ^8f by solving the BS 
equation with a realistic one-boson exchange kernel including tt, u, p , a, rj, <5 mesons. In 
the adopted approximation, the BS amplitude reads (see also Appendix |AJ) 


* M {k) = *&+(k) + *%++&). 


M 


(16) 


The charge-exchange vertex A ce is also a matrix in the spinor space and, for consistency 
of the approach, it would be preferable to decompose it into partial vertices as in the case 
of the BS amplitude and to find the coefficients from the NN charge-exchange reactions. 
However, this is a rather cumbersome procedure in which one can determine only the 
vertices for on-shell NN processes. To extend it to off-shell nucleons one needs to specify 
some method, which inevitably requires theoretical models and additional approximations. 
In our case, in the process O the transferred momenta and energies are considered 
relatively small, hence the virtuality of nucleons in the vertex A ce may be disregarded. 
Then A ce may be expressed directly via the amplitude f r 's',sr of the real charge-exchange 
processes p + p n — p p + n with all particles on the mass shells 

fr's',sr = ui ( p p ) U r p (n) A^g U^{p) u|(p n ). (17) 


Then 


T m = y- 

r'r / j 


1 (2m) 


d A k f r 


x u s (jp n )^ M (k) (^D -k + m)^ P '(k - ^q) u s \p p ). 


(18) 


Note that the amplitude (TH) is manifestly covariant. For the final 1 S'o-state within the 
p-spin classification the BS amplitude in the center of mass of the NN pair is represented 
by four partial amplitudes 1 S' ( )" + , 3 Pq~ and 3 P 0 ~ + |I4|, which for the sake of brevity 

in what follows are denoted as 0i,...,0 4 . In order to avoid an explicit Lorenz boost 
transformation to the laboratory system, it is convenient to write the 1 S'o amplitude in a 
covariant form 

1 


a/47T T pi (p) = 

- 63(75 


6175 - b 

pi — m 


’2 — (75Pl +P275) 
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P2 + m . , p2 + m pi—m 

-75) - 64—7^—75- 
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(19) 
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where p 12 = P'/2 ± p , and p is the relative momentum. The four Lorenz invariant 
functions fe* = bi(P'p,p 2 ) in the center of mass of the pair are linear combinations of the 
amplitudes (pi = (pi(r 0 , |r|), i = 1, ...,4 j34| (see Appendix |B[). Now it is sufficient to 
express the amplitude of the process ([[]) in terms of deuteron ” + +” components and (|l^) 
to implicitly account for the Lorenz boost effects [IT]. For the final state of the pair, as 
in the deuteron case, all the amplitudes with negative p-spins are neglected as well. 

Substituting eqs. © and ([T6|) into eq. flTE|) the matrix element may be written in 
terms of two-component spinors and 3-vectors as 


u s {p n ) d 'M(k) k + m\ ^!p, [k - 0 u s \p p ) 


m \M d — k 0 — E 


' £m) Xs 


'*-%h 


16tt E e 

+ m ) ( Pp + ?«) 

+xl{*r ' k ) Xs'( k • €m)^dC 2 

+ -xl(*r ■ q) X,' ( k • £m) (V s + V2*p D ) + xl(c t ■ k) Xa> (q ■ £m) 

ip d - 


+xlxs'i([<l X k] • £a/) tys 




c 3 > M r o, |r|)), 


( 20 ) 


where the quantities Ci, C 2 , C 3 have a pure kinematical origin and are independent of spin 
variables. Their explicit form can be found in the Appendix |E| Now, from eq. (p0|) it is 
clearly seen how to compute Epl at given spin variables r', r, M and, consequently, the 
invariant amplitudes (|ll|) and observables (|12|) , ([131). The partial amplitudes 

(pi may be found from the BS equation, which, in the simplest case of pseudo-scalar 
exchanges reads as 

d A p' 


T P ,(p) = T p ,(jj) + ig^ NN 


(2tt) 


■ A(p - p')S(p 2 )^p'pp'^Sppx), 


( 21 ) 


where A and S are the scalar and spinor propagators respectively, S = UcSUq 1 , and 
Tp, (p) is the free amplitude corresponding to two non-interacting nucleons (the relativistic 
plane wave). The solution of eq. (^l|) may be presented as a Neumann-like series, the 


first term of which is the first one from eq. (21): 


tp.(p) = 


( 22 ) 


The second part in eq. (|22|) is entirely determined by the interaction and may be sym- 
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bolically referred to as scattered wave. For the 1 5' 0 state one has 

y°Ar)\p'= {V sj fi) = |r|)ri 5 ++(r), 


0?(>o, M) = 2(2t t) 4 


1 1 


5(r 0 )5(|r| - |r*|), 


y/An | r* | 2 

where r = (ro, r) is the relative 4-momentum of the pair (the variable of the momentum 
space), |r*| = ^Sf/A — m 2 is the experimentally measured relative 3-momentum of the 
is the spin-angular harmonic for the 1 S'o state [|54 


pair, 


is++W 


To determine the 


scattered wave in eq. (|2|) it is necessary to solve the BS equation of the type ([Hp 
including all the above mentioned exchange mesons. Solving the BS equation in the 
continuum is a much more cumbersome procedure than solving the homogeneous BS 
equation. Besides difficulties encountered in solving the latter (singularities of amplitudes, 
poles in propagators, cuts etc.) the former even does not allow the usual Wick rotation 
|3?H to the Euclidean space, and there are no rigorous mathematical methods to solve eq. 


( PH) in the Minkowsky space |48|]. However, an approximate solution of eq. (|2l|) may be 
obtained by employing the so-called ’’one-iteration approximation” |TT| , within which one 
may obtain a rather good estimate of the interaction term (see below). 


IV. RELATIVISTIC IMPULSE APPROXIMATION 

We start our analysis of the reaction (|l|) by disregarding the interaction term in eq. 
0, i e. putting 

0i(r o , |r|) = 4>i(r 0 , |r|), (p 2 = 0, 0 3 = 0, (j) 4 = 0. (23) 


In this case the final state of the pp pair is described by the free part 0?(r o , |r|), what 
obviously means the 1 S'o part of two plane waves. Within the impulse approximation all 
formulae become much simpler and one may preliminary investigate the main features of 
the process ([!]), fix parameterizations of the elementary charge-exchange amplitude, find 
proper kinematical regions where the assumptions made hold etc. 

One has: 


S(n) = 6 l(P',k--q)/Sp], 


(24) 


<5(|r| — |r*|) = 5 


ik--q ) 2 


i s f 


m* 


(25) 
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and 


|r|) = 2(2t t) 4 


s f 


^ .£7|k| |q| V \sf - m 2 


x ,)ii Vl - [1 m d - £]] i (cos 0j :<l + i 2 ^ ( k ('(° * ) 


(26) 


In eq. (]2 (d|), 9k q is the angle between k and q. Then, with eqs. (|2B|), ( pj| ) the matrix 
element in (1X81) reads 


max 

/ , 

J Md 

J 

V 27T |q| 


|k|d|k| 


2t r 


d(pk fr's',s 


X |xi(o- • £m)Xs' (u s - Qs/ + mP' Q 

+ xi(<r ■ RqK)Xs'{£, M ■ R q k)A=U D P ° 


y/2 E — m 

- xli* • q)Xs'(^M • R qk) (u s + y/2U D ) + ■ Rqk. )Xs'(q ■ £m) (Vs - 


+ <S M »i([q x P g k] ■ £ M ) ^7 S - 


tV 


(27) 


(28) 


Here we introduced the corresponding BS wave functions 

7 j = G S,P _ 

S ’ D 47T \/2M d {M d — 2E) ’ 

where are the BS vertices pjfl . In this notation, the introduced quantities Us.n 

correspond fully with the non-relativistic S and D wave functions of the deuteron, and 
the non-relativistic treatment of the results in terms of usual wave functions becomes 
more transparent. In eq. (^), the matrix R q describes the rotation about the y axis by 
an angle 6 (the angle between q and z axis) defined as 


Rq = 


^ cos 9 0 sin 9 \ 
0 1 0 
\ —sin 9 0 cos# / 

and the limits of integration over |k| are as follows 


max, min 


N 


i + 


q- 

s f 


\sf-m 2 ± ||q| 
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Now it is straightforward to compute T 7 y/ for any value of spin indices and, by virtue 
with eq. (0), the quantities M, £?,...JF. Here it is worth commenting how the charge- 
exchange amplitude f r ' s ',sr is involved into our numerical calculations. Beside spin indices 
this amplitude also depends upon two Mandelstam invariants being the total energy of 
the subprocess of charge-exchange, s pn = (D/2 + k + p) 2 (where p is the momentum 
of the initial proton), and the invariant transferred momentum t = (n — p ) 2 , which is 
a common variable for both the full reaction (il|) and the subprocess of NN reaction. 
From kinematics, the 4-momenta of the ’’initial” neutron and ’’final” proton are D/2 + k 
and P '/2 + k — q/2, respectively, which might be both off-mass shell. In the relativistic 
impulse approximation after integration with the 5 function, the final proton receives an 
on-mass shell momentum, p p = P '/2 + k — q/2. Consequently here only the neutron from 
the deuteron remains off-mass shell. Hence we are left with a charge-exchange amplitude 
with one nucleon off-mass shell and a varying yfs. In ref. |29|, the dependence of the 
charge-exchange amplitude upon the initial energy has been found to be rather weak 
(at initial energies in a range of few GeV). Therefore in our calculations, by neglecting 
the off-mass shellness of the neutron, the charge-exchange amplitude is taken from the 
real NN process at equivalent values of yfs (the invariant energy of the incoming proton 
and internal neutron). Moreover, since this amplitude is essentially independent of the 
azimuthal angle it is taken out from the corresponding integration. Such a procedure 
of implementation of real amplitudes into calculations where one or several nucleons are 
off the mass shell is commonly adopted in literature |14], |39|, [40], [41| with an ” a posteriori” 
justification from the comparison of numerical results with experimental data. 


A. Numerical results 


In our numerical calculations we employ a parameterization of the elementary charge- 


exchange amplitude from ref. [^9| and parameterizations from partial-wave analysis per¬ 
formed by different groups 


1143[], which are available via web, cf. 1 |44| . [45| | . The partial 
charge-exchange amplitudes are available in the helicity basis as partial helicity ampli¬ 
tudes 


/ 1 = ( + + |M ce | + + ), / 2 = ( + + |M ce |--), 

/ 3 = ( + - |Fl ce | + -), 
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f 4 = ( + - \A ce \-+), h = ( + +\A ce \ + -), 


(29) 


normalized as 


do ct 

dt 


327rs(s —4m 2 ) \~[ 


X l/*l 2 + 4|/ 5 


(30) 


Since in our matrix element ( |27|) the spin amplitudes f r ' s ',sr are defined in the deuteron 


center of mass, the hclicity amplitudes (|29|) must be first boosted along the direction p + k 
from the center of mass of the pn system to the laboratory system, and then transformed 
to spin amplitudes by Wigner rotations. Taking into account that the procedure of Lorenz 
boost itself for each particle results in an additional hclicity Wick rotation ||46]], one needs 
eight rotations for each amplitude in eq. (|29|) (see Appendix 0 for details). 

In the Fig. [2] the partial hclicity amplitudes ([0) are presented as a function of the 
transferred momentum |q| and energy corresponding to the initial momentum in the lab¬ 
oratory system |p| = 2.5 GeV/c. The full lines represent the partial-wave analysis ||3 


whereas the dashed line depict results of the analytical parameterization from ref. 

Both parameterizations describe equally well the unpolarized charge-exchange cross sec¬ 
tion, nevertheless a substantial difference is seen between two sets of parameterizations. 
Obviously, the unpolarized cross section is not sensitive to details of partial amplitudes, 
and for a more precise determination one needs more independent measurements of po¬ 
larization observables. In this context, we remark that such reactions can be considered 
as an additional source of information about the elementary amplitude, as pointed out, 
e.g., in ref. [|3|], since the process ([!]) is entirely governed by the elementary subprocess 


of NN charge-exchange. 

In what follows we are interested in systematical calculations of the cross section 
and tensor analyzing power of the process (jlj) for kinematical conditions achievable at 
COSY [Q. However, first we perform calculations for such kinematical conditions for 


which experimental data are already available [|19|]. Note that experimentally one mea¬ 
sures the cross section averaged over some interval of the excitation energy of the pp pair. 
Conformly, we define the differential cross section do/dt as the double differential cross 
section (|6j) averaged over a given bin of energy 


do' 

dt 


(87t) 3 A 


dsf 


Rk 


\ 


1- —IM/.I 2 . 

S f 


k = 1,2,3..., 


(31) 
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where k labels the intervals of E x given in the experiment. At SATURN-II ]19|, where 


the process (|1) has been investigated in details at initial momenta of protons p 
and 0.599 GeV/c, the mentioned intervals of E x are 

= 0.444 

R\ : 

0 < E x < 1 MeV, 

(32) 

R 2 : 

1 < E x < 4 MeV, 

(33) 

R 3 : 

4 < E x < 8 MeV. 

(34) 


The intervals R\ and R 2 fit into the COSY kinematics H as well. Note, that under the 


kinematical conditions (|32|-[f|) the variable t is indeed small, ranging in an interval from 
0 to 0.16 (GeV/c) 2 . 

In Figs. [| and [| the cross section and tensor analyzing power T 20 for the process © 
are presented. The initial energy corresponds to a typical COSY momentum |p| = 2.5 
GeV/c. The solid lines depict results with the elementary amplitude taken from ref. [fh| 


whereas the dashed lines are results with the parameterization from ref. |29|. As expected, 
since different parameterizations equally well reproduce the elementary charge-exchange 
cross section, the unpolarized cross section of the process (jT|) is not sensitive to parame¬ 
terizations of the elementary amplitude. An opposite situation occurs when calculating 


the polarization observables defined in eq. (|T2|) for which the contribution of partial am¬ 
plitudes is non-diagonal and interferences might be important. This is clearly illustrated 
in Fig. [| where the tensor analyzing power (|T3| ) exhibits indeed a strong sensitivity to 
parameterizations of partial amplitudes and practically does not depend upon the chosen 
bin of the excitation energy. From this picture one may conclude that an experimental 
investigation of the tensor analyzing power may constrain further parameterizations of 
the elementary charge-exchange amplitude at high energies. 

As already mentioned, the process © has been experimentally investigated at 
SATURN-II 0. Although at such energies the final state interaction cannot be ne¬ 
glected and the simple impulse approximation is too rough, a comparison of data with 
theoretical results is rather instructive. In Figs. and |6] we present results of calculations 
of the cross section (0) and tensor analyzing power T 2 0 ( |T3| ) defined by eqs. (|27|-pl|) 
together with available experimental data. The full lines have been obtained with param¬ 
eterizations of the elementary amplitude from refs. [|3|, [I3|, while the dashed lines with 
the parameterization given in ref. [^9|. From Fig. it is seen that the impulse approxi- 
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mation qualitatively describes the general shape of the cross section as a function of the 
transferred momentum |q|. For low values of |q|, say up to 0.2 GeV/c, and in the interval 
of the pair excitation energy 1 < E x < 4 MeV here is even a good agreement with data, 
in contrast with other intervals of E x and higher values of |q|. From this and from the 
results of non-relativistic calculations [H], where final state interaction and higher partial 
waves have been taken into account, one may conclude that at high values of transferred 
momentum the effects of final state interaction in the 1 S'o state become dominant. At 
higher excitation energies the interaction effects are not so significant, however here cor¬ 
rections from other partial waves may become important. The same conclusions can be 
drawn from Fig. ]Bj. where the tensor analyzing power, computed with two parameteriza- 
tions (as above, the solid lines correspond to ref. (42], [43]], dashed curves to ref. S). 
is compared with experimental data. From Fig. [| it is also obvious that a qualitative 
agreement with data for the tensor analyzing power may be achieved only by using the 
elementary charge-exchange amplitude from the partial-wave analysis (42|, [43[], while the 
parameterization |29|] results even in opposite sign for ^20- This is a direct indication that 
a more sophisticated partial-wave analysis gives more reliable partial helicity amplitudes. 
Nevertheless, since such an analysis has been performed for low and intermediate energies 
(up to few GeV), a further tuning of partial amplitudes (|29|) at relativistic energies is 
still desirable. Together with the proper choice of kinematics at high energies (i.e. a 
kinematical situation where the role of higher partial waves in the final state, e.g., triplet 
states, is suppressed |23|]) one may expect that the proposed mechanism will adequately 
describe reactions of the type (|T;) and corresponding information may be obtained. Note 
that in all the above calculations the vector analyzing power is strictly zero. 


V. ONE-ITERATION APPROXIMATION 

As mentioned, for a consistent relativistic analysis of reactions with deuterons and two 
interacting nucleons in the continuum one should solve the BS equation for both bound 
state and scattering state within the same interaction kernel. We have found a numerical 
solution for the deuteron bound state with a realistic one-boson exchange potential |(27j1 . 
The BS equation, after a partial decomposition over a complete set of matrices in the 
spinor space, has been solved numerically by using an iteration method. We found that 
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the iteration procedure converges rather quickly if the trial function is properly chosen, 


e.g., if in the BS equation the combination of the type (|28|) is used as trial functions 
with non-relativistic solutions of the Schrodinger equation. In such a case, even after the 
first iteration, the BS solution coincides with the exact one up to relative momentum 
p ~ 0.6 — 0.7 GeV/c. This circumstance can be used if one needs an approximate solution 
of the BS equation at not too large momenta p < 0.5 — 0.7 GeV/c. This is just our 
case, since in reaction (|T|) the relative momentum of the pp pair is expected to be rather 
small and the scattering part of the amplitude (|22|) can be obtained from eq. ( 0 ) by one 
iteration, provided the trial function is properly chosen. 


Formalities 


To solve eq. (0) we proceed as follow (cf. ref. [in]]): (i) for simplicity, in the inhomo¬ 
geneous BS equation we leave only the pseudo-scalar isovector exchange (tt- mesons), (ii) 
write instead of eq. fl2ll) the mixed BS equation by introducing in both the left hand side 
and the free term the BS vertices, i.e. 

d 4 p' 75 v hp(p / )75 


Gp{p) — G° P (p) J 


(35) 


(2v i) 4 {p - p'Y - gV 

(iii) bearing in mind that the BS partial vertices may be obtained from the same spin 


angular functions r a (p), by replacing p —p |36|, we write the corresponding partial 
BS equation 


GV+(po, |p|) = G° ls ++(po, |p|) 


. 2 f ddp'dVtp (E p E p , - m 2 + (pp')) 0i(p' o , |p'|) 

l 9irNN 


(36) 


- (27T)- E r Ey (r-r/y-tf, 

Further by disregarding the dependence upon po in the meson propagator in eq. ( |36| ) and 
then using the standard representation of propagators via generalized Legendre polyno¬ 
mials Qi and restoring the BS equation in terms of partial amplitudes, one obtains 

1 


9nNN 


* sr (P0, IPI) = 0V(PO, IpI) - ^77- v 

V2 v*/ ^p) 


X 


Po 


°? 71 / 

f dp 


P' 

1 r 

J 2 tt 

0 

P 



(■ EpEp' - m 2 )Q 0 (yp) - |p| \p\Qi(y„j\ ui So (s f , |p'|), (37) 
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p 2 p 2 + /i 2 

where y M = -——-— -. In obtaining ([57]) the integration over p' 0 has been carried out 

21P11 p' | 

1 


in the pole p 0 = -y's/ ~ E p' an< ^’ si m il ar to eq. (|28|) , we define the BS wave function in 
the continuum as 


ui So (s f , |p'|) = 


9i S ++(Po, Ip'I) 
V^f- 2E P' ' 


(38) 


Now, if we restrict ourselves to only one iteration in (|37|) taking the trial function (f38|) 
as a non-relativistic solution of the Schrodinger equation, e.g. the Paris wave function 
u[ s^( s /’ Ip'I), th e BS amplitude is obtained as 

G oA \p o, |p|) 


015++(Po, |p|) = <f>is++(Po, IpI) - . N 2 

-pS 

where the ’’one-iteration” BS vertex G OA '(po, |p|) is defined by 


(39) 


G oi (po, |p|) = 

^T 47T 


E v 

1 -^ 

m . 


dr e ^jo(pr)uu;f(r) 


iPl 


mE„ 


• u NR (r) 

dr _±I2 —L e ~P r (i _ ^ r ) jd(p'r) } . 


(40) 


From eqs. (39) and ([h]) one can easily find the non-relativistic analogue of the obtained 
formulae. The free term in eq. (^) together with the first term in eq. (|40|) reflect the 
non-relativistic equation for the 1 S'o wave function, while the second term in ([ID]) turns 
out to be a correction of purely relativistic origin. 


B. Numerical results 


In Figs. 0 and we present results of numerical calculations of the cross section and 
tensor analyzing power T 2 o given by eqs. ([31]), (|T3D, (H), ( [20] ) and (p9|-|40|) . The elementary 
charge-exchange amplitude has been taken from ref. |43| and the non-relativistic trial 
function u^(r) as the solution of the Schrodinger equation within the Paris potential [[47] . 
The BS S ++ and D ++ amplitudes are those from the numerical solution |27| obtained with 


a realistic one-boson exchange interaction. The dashed curves in Figs. [7j and [R| correspond 
to results within the relativistic impulse approximation, while the solid lines depict results 
with taking into account the final state interaction in one-iteration approximation. It is 
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seen that in all three energy bins the agreement with data for the cross section is essentially 
improved. This concerns especially the range 1 < E x < 4 MeV. For the energy bin close 
to zero there is still a disagreement with data at low transferred momenta which probably 
may be related to the fact that in onr calculations we have not taken into account the 
Coulomb interaction within the pp pair. For higher excitation energies (E x ~ 8 MeV), 
other partial waves (e.g. triplet state) in the pp final state contribute and, within the 
adopted assumptions, one may expect only semi-quantitative agreement with data. From 
Fig. 0 one may conclude that at low excitation energies the supposed mechanism for the 
reaction (HP (i-e. charge-exchange subprocess with interaction in 1 S 0 state of the pp pair 
in the continuum) seems to be correct. Moreover, from a comparison of the left and right 
panels in Fig. [0 one may expect that the higher initial energy the larger kinematical region 
where the mechanism holds. Fig. |] demonstrates that the tensor analyzing power is less 
sensitive to final state interaction effects. As a matter of fact, the tensor analyzing power 
(|T3|), being a ratio of non-diagonal products of partial amplitudes to the diagonal ones, 
serves as a measure of the quality of parameterization of partial amplitudes and their 
mutual relative phases. This has been pointed out in a series of publications (see e.g. 
refs. Ii), where a good simultaneous description of cross sections and T 20 in reactions 
of the deuteron break-up or elastic scattering of protons, is still lacking. Nevertheless, 
since in the process (P the behavior of the partial amplitudes (|TT|) , as seen from eq. 
(POj), is mostly governed by the elementary charge-exchange amplitudes, an experimental 
investigation of the tensor analyzing power T 2 q in reactions of the type (P can essentially 
supplement data on the NN charge-exchange amplitudes at high energies. 

In Figs. H and |T(] we present the predicted cross section and tensor analyzing power at 
high energies relevant for COSY and Dubna accelerator. It is immediately seen that the 
cross section is substantially decreasing with increasing energy, nevertheless it remains 
large enough to be experimentally easily accessible. Another peculiarity of the studied 
process at relativistic energies is that the tensor analyzing power T 20 does not change the 
sign remaining positive in a large kinematical region, in contrast to lower energies (cf. 
Fig. H). Note again, that in the above calculations the vector polarization of the deuteron 
is strictly zero. 

From the performed analysis one can conclude that there is a kinematical region for 
the excitation energy, E x < 5 MeV, and transferred momentum, |q| < 0.3 — 0.4 GeV/c 
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(i.e. the COSY |Q kinematics), for which the mechanism of the reaction (|lj) is fairly well 
described within the spectator approach by an elementary pn charge-exchange subprocess, 
for active nucleons, with detection of the pp pair in the 1 S'o final state. Our covariant 
approach agrees with previous non-relativistic calculations and allows for predictions of 
the cross sections and polarization observables at intermediate and relativistic energies, 
in particular, for kinematical conditions which are realized at COSY. The predicted cross 
sections of the process and the tensor analyzing power T 20 are large enough to be used, 
in a large range of initial energies, for determining properties of the polarized deuteron, 
provided experimentally one simultaneously detects a vanishing vector polarization of 
deuterons. 


VI. SUMMARY 

In summary, the performed covariant analysis of the reaction D(p,n)pp with the two 
final protons in a 1 S'o state allows us to conclude that, as in the non-relativistic limit, such 
a process can be used as an effective deuteron polarimeter also at relativistic energies, 
in particular, at the range covered by COSY at Jiilich and upgraded Dubna accelerator. 
Additional information about the elementary charge-exchange amplitude at high energies 
can be obtained from precision data with known deuteron polarization. 
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APPENDIX A: DEUTERON STATE 


The BS amplitudes of the deuteron, eq. ([16]) , in the deuteron center of mass explicitly 
read 


^f++(k) = A/'(fci + m ) 1 + 2 70 i M (k 2 - m)ip s (ko, ]k|), 


(Al) 


a r 


1 + 7o 


®D++( k ) = 2 

3 


(A 2 ) 


x (h - k 2 ){k^ M )j (h ~ m)^ D (k 0 , |k|), 

where k\ t2 are on-mass shell 4-vectors, 

k\ — (E, k), k 2 — (E, -k), k — (k 0 , k), E = \j\c + m 2 , (A3) 

and ips,D(k 0 , |k|) are the partial scalar amplitudes, related to the corresponding partial 
vertices as 


fpsMko, |k|) = 


Gs,D(k 0 , |k|) 


M, 


D 


E) -kl 


In Eqs. (|A1[-|A2|) the normalization factor is A f = 


^/8tt2E(E + m) 


APPENDIX B: 1 S 0 STATE 


Different representations of the BS amplitude 1 S'o in the continuum have been studied 


in details in ref. [34]], where the reader may found the most general expressions for the 
covariant amplitudes b t , eq. fll9l) , in terms of the partial amplitudes (j)j in the center of 
mass of the pp pair. Since in the present paper we consider only the ++ component of 
the (j)i , we are left with one invariant function b , which is taken to be 64 . Then the explicit 


expressions for the kinematical coefficients C±,C 2 ,Cs in eq. (| 20 | ) can be cast in the form 


Ci = (p ■ k) Ai - (E + m) ( p° + m) A 2 - (q ■ p ) (E + m) A 3 - ( q • k) (p° + m) A 4 , 


C 2 = 


\/2(E- 


m 


[ ( E-m)A 1 - (p° + m) A 2 — (( 2 k - q) ■ q) A 3 ], 


C 3 = Ai - (E + m ) A 3 + (pp + m ) A 4 , 
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where 


M = K, [Pq (p° p + 2m) + (q • p p )] - 1, 
M = K [Pq ( p° p - 2 m) + (q ■ p p )] - 1, 
M = -/C(2m-p° + pJ), 

A 4 = -K,(2m + p° 1 -p 0 p ), 


(Bl) 

(B2) 

(B3) 

(B4) 


with the invariant coefficient /C 
by 


(see ref. ]04| ). The 4-vectors p\ and p p are defined 


Pi = \ p ' + k ~ \q= (P?,P P ), 

Pp = (pj, P p ), P° p = \J r n 2 + Pp, Pp = k - q. 


Observe, that p p represents an on-mass shell vector. Note also, that in the relativistic 
impulse approximation, since p\ = p p , eqs. ( [BT| - [B4l) are substantially simplified. 


APPENDIX C: RELATIVISTIC SPIN TRANSFORMATIONS 

By definition, a state with given momentum p and hclicity A in a frame of reference O 
is that obtained by a Lore n z transformation of a state with given spin projection s z from 
the rest system O rest to O , i.e.: 


|p; A) = | P,s,8 g ) 0 , (Cl) 

where P— (m, 0, 0, 0). As usual, a Lorenz transformation h[ p] is presented by a sequence 
of two operations: a boost along the z axis, l z (v), where v is the speed of the state in O, 
and a rotation from z direction to the direction of p, i.e. O = r^ x (0, 9 , 0 )l~ l {y)O rest . 

Let us suppose now that one has a state |p; A) given in the frame O and one wishes to 
know how it reads in another frame O' obtained by a Lorenz transformation l on O 

|p ■X) ey = U(l- 1 )\p;X). (C2) 

From the definition of the helicity states one has 

C/(r')|p; A)=C/(r 1 )C/(/i[(p)])|P;A), 


(C3) 
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where h[p] is the corresponding Lorenz transformation P —> p. Then multiplying eq. (|C3l) 
by unity, U(h[(p')])U~ l [h[(p')] = 1, where h[ p'] is the hclicity transformation that dehnes 
a state |p'; A) = f/(h[(p'])| P, A) with p' being the same vector as obtained by transforming 
p from O to O', one obtains: 

u(i- 1 )\p:\) = u(h{( P f ])n\p,\), (C4) 

where 1Z = f/ _ 1 [h[(p , )]f/(/ _ 1 )f/[h(p)] is the sequence of transformations P—> p —> p' — 
i.e. nothing but a rotation. Then, 

|p,A> 0 . = B«H|p',A'>, (C5) 

where u is a set of Euler angles describing the rotation. In case when the Lorenz transfor¬ 
mation is a simple boost along the z direction with the speed /3, then u is just an angle, 
describing a rotation about the Y axis, 

cos u = cos O' cos 0 + 7 sin O' sin 6 , (C6) 

with 7 = l/\/l — /3 2 , and 0,0' are the polar angles of p in the systems O and O', re¬ 
spectively. This is known as Wick helicity rotation, contrary to Wigner’s canonical spin 
rotation. In our case, the relevant z axis is the one along the direction of (k + p). Then, 
obtaining the helicity amplitudes in the laboratory frame we need an additional rotation 
to change from the hclicity basis to the spin projections. 
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P 11 


FIG. 1: Spectator mechanism for the charge-exchange process pD —> n(pp). The Bethe-Salpeter 
amplitude for the deuteron bound state and the pp pair in the continuum are denoted as and 
T, respectively. The elementary pn charge-exchange amplitude is symbolically represented by 
„ 4 ce . 
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FIG. 2: Partial helicity amplitudes of eq. (|2fj) vs. the transferred momentum |q| for different 
parameterizations. Solid curves correspond to the partial wave analysis of |43|, 45], while the 
dashed curves are the high energy parameterization given in [p^]. The amplitudes are dimen¬ 


sionless. 
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|p| = 2.5 GeV/c 



FIG. 3: Unpolarized cross section, eq. (31), calculated within the relativistic impulse approxi¬ 
mation at a typical COSY initial momentum for several bins of excitation energy of the pp pair. 
Results of calculations with two parameterizations of the NN charge-exchange amplitude are 
exhibited (solid curves are obtained with the amplitude from |43|, [hj, the dashed curves from 
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|p| = 2.5 GeV/c 



FIG. 4: The same as in Fig. ^ but for the tensor analyzing power T 20 (13). 
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|p| = 444 MeV/c |p| = 599 MeV/c 



|q|, GeV/c 


FIG. 5: Comparison of the non-polarized differential cross section eq. (|3l]) within the impulse 
approximation at non-relativistic initial energies (|p| = 0.444 GeV/c (left panel) and |p| = 0.599 
GeV/c (right panel)), with experimental data [|1|. Solid curves correspond to the elementary 
charge-exchange amplitude from the partial wave analysis from refs. [^3|, ^], dashed curves use 

Hi- 
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|p| = 444 MeV/c |p| = 599 MeV/c 



|q|, GeV/c 


FIG. 6: The same as in Fig. |] but for the tensor analyzing power T 20 (|l3|). 
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|p| = 444 MeV/c |p| = 599 MeV/c 



|q|, GeV/c 


FIG. 7: Results of calculations of the differential cross section eq. (|31|) with taking into 
account the effects of final state interaction in 1 So state (solid curves). Experimental data are 
those from SATURN-II [19], the elementary amplitude has been taken from refs. J43|, ^5|. The 
dashed curves depict the results of calculations within the pure impulse approximation (cf. solid 
curve in Fig. 5). 
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|p| = 444 MeV/c |p| = 599 MeV/c 



FIG. 8: The same as in Fig. |?] but for the tensor analyzing power (13). 
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|p| = 2.68 GeV/c 



|q|, GeV/c 


FIG. 9: Results of calculations of the differential cross section ( |3l| ) with taking into account the 
effects of final state interaction in 1 Sq state (solid curves). Kinematical conditions correspond to 
those proposed in Q for experiments at COSY. Elementary amplitude from refs. ^5|. The 
dashed curves depict results of calculations within the pure impulse approximation (cf. solid 
curves in Fig®. 
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|p| = 2.68 GeV/c 



FIG. 10: The same as in Fig. but for the tensor analyzing power (13). 




